Activated bond-breaking processes preempt the observation of a sharp glass-glass 
transition in dense short-ranged attractive colloids 
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We study — using molecular dynamics simulations — the temperature dependence of the dy- 
namics in a dense short-ranged attractive colloidal glass to find evidence of the kinetic glass-glass 
transition predicted by the ideal Mode Coupling Theory. According to the theory, the two distinct 
glasses are stabilized one by excluded volume and the other by short-ranged attractive interactions. 
By studying the density autocorrelation functions, we discover that the short-ranged attractive glass 
is unstable. Indeed, activated bond-breaking processes slowly convert the attractive glass into the 
hard-sphere one, preempting the observation of a sharp glass-glass transition. 

PACS numbers: 64.70.Pf, 82.70.Dd 
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Short-ranged attractive colloidal systems have re- 
cently become the focus of many experimental 1] and 
theoretically studies. The interest in these systems stems 
from their peculiar dynamics for showing structural 
arrest phenomena both of gelation and vitrification type, 
and for being amenable of analytic treatments. Previous 
studies have convincingly shown that unusual dynam- 
ical phenomena emerge from the competition between 
two characteristic localization length scales; the hard- 
core and the short-ranged attraction localization lengths. 
In hard-sphere colloids, when the packing fraction (f) ex- 
ceeds roughly 0.58, the colloidal suspension turns into a 
glass and particles are confined in cages of size of the 
order of 10% of the hard core diameter (the hard-sphere 
localization length). In attractive colloids, localization 
can be realized not only via excluded volume (i.e. tuning 
(/)) but also thermally (i.e. tuning the ratio between tem- 
perature T and potential depth) Hence, slowing down 
of the dynamics can be induced not only by increasing 
but also by progressively lowering T, i.e. increasing the 
inter-particle bonding. Previous studies have shown that 
when the range of the attractive interaction is smaller 
than 10% of the hard-core diameter, an additional local- 
ization length, controlled by the potential, sets in and 
an efficient mechanism of competition between the hard- 
sphere and the short-ranged bonding localization lengths 
arises. In the fluid phase-region where the two localiza- 
tion mechanisms compete, a highly non trivial dynam- 
ics is observed: particle diffusion shows a maximum on 
heating, the fluid can be transformed at constant </) into a 
glass both on cooling and heating, the time dependence of 
dynamical quantities like the mean squared displacement 
and the density-density correlation function, have an un- 
common behavior at intermediate times, showing respec- 
tively a sub-diffusive and a logarithmic regime. These 
anomalous properties have been recently observed in a 
series of beautiful experiments [1, IE 13 > and exten- 
sive simulations of particles interacting via short-ranged 
so uare- well jili 1 1 Ol 1 1 If and Asakura-Oosawa potential [l3|. 

Hard sphere colloids have been an important system 



for accurately testing theoretical predictions. It has been 
shown that the ideal Mode Couphng Theory (MCT), 
which neglects hopping phenomena, provides an accu- 
rate description of the dynamics close to the hard-sphere 
glass transition^E^3|- Ideal MCT provides an accurate 
description of the fluid phase dynamics also in short- 
ranged attractive colloids, despite the complex dynam- 
ical processes alluded to previously. Several theoretical 
predictions [H, HE 113, 13 E 113 have been recently 
quantitatively confirmed by experiments |E IE 13 Hi and 
simulations ,6^.9. .10llllLlia|. 

An untested important prediction of the theory regards 
the presence of a kinetic (as opposed to thermodynam- 
ics) glass-glass transition which should take place in the 
glass phase on crossing a critical temperature. Heating 
a short-ranged attractive glass should produce a sudden 
variation of all dynamical features, without significant 
structural changes. For example, at the transition tem- 
perature, the value of the long-time limit of the density- 
density correlation function, the non-ergodicity factor fq, 
should jump from the value characteristic of the short- 
ranged attractive glass to the significantly smaller value 
characteristic of the hard-sphere glass. 

In this Letter we report an extensive numerical study 
of the binary short-ranged square well (SW) potential in 
the glass phase aiming at studying the T-evolution of the 
glass dynamics and detecting the location of the glass- 
glass transition. We discover that hopping phenomena 
destabilize the attractive glass and prevent the observa- 
tion of the a discontinuous glass-glass transition. 

We perform molecular dynamics (MD) calculations of 
a 50% — 50% binary mixture of = 700 particles of 
unitary mass interacting via a square-well potential of 
unitary depth —uq. The asymmetry between the two 
hard-core diameters is fixed to 20%. The attractive well- 
width Aij is given by Aij/(Ay + aij) = 0.03, where Cp- 
is the hard-core diameter for the ij-type interaction|2l|. 
with i,j — A,B. In this model, a bond between two 
particles can be unambiguously defined when the pair 
interaction energy is —uq- The one-component version 



of this model has been extensively studied within MCT 
|l7|. The binary mixture case, which allows to prevent 
crystallization without substantial changes in the dynam- 
ics, has been studied numerically in the fluid phase [loj . 
The iso-diffusivity curves in the cj) — T region are re- 
entrant and, close to the re-entrance, correlation func- 
tions show a logarithmic decay in agreement with MCT 
predictions. Comparing the theoretical MCT calcula- 
tions for this specific model, with the numerical results 
reported in Ref. ^3 — following the procedure first used 
by Sperlji^l — the glass-glass transition line has been lo- 
cated between {(f) ~ 0.625, T ~ 0.37) and the end-point 
{(t> ~ 0.64, T ~ 0.42) Pi . At = 0.635, the packing 
fraction studied in this work, the glass-glass temperature 
is predicted to be T ~ 0.4. 

To generate glass configurations, we compress 24J, to 
4> = 0.635, 136 independent equilibrium configurations 
generated at T = 0.6 and (j) ~ 0.612 in a previous 
studv[lo|. We checked that the results presented in this 
Letter do not depend on the T and (f> of the starting un- 
compressed configurations. Since the simulated system 
is out-of-equilibrium, ensemble averages are requested to 
reduce the noise level. Each compressed configuration 
was instantly quenched to several temperatures T by 
proper rescaling of the initial velocities and let to evolve 
in a NVT simulation |2^J. On the time scale of our study 
the system behaves as a glass at all T values, i.e. no diffu- 
sional processes are observed. We analyze nine different 
T values, from T = 1.5 to T = 0.1. Each T requested 
about 15 CPU days on fast AMD Athlon processor. 

Since the system is in an out-of-equilibrium state, cor- 
relation functions depend not only on t but also on the 
time elapsed from the compression, the waiting time t^. 
For t « tu,, correlation functions for different tw col- 
lapse on the same curve — a phenomenon connected 
to the equilibration of the degrees of freedom faster 
than iu;[2a|. We exploit such feature to generate t^- 
independent data and focus only on the T-dependence. 
For this reason we consider the dynamics for times up to 
approximately 400 time units and > 4000. 

To confirm that in the time window < t < 400, no 

-dependence is observed for > 4000, we calculate 
the density-density correlation function 




FIG. 1: Waiting-time dependence of (jtqi'tw + t,tw) for T — 
0.35, at qoBB — 25. For tw > 4000, no aging is observed in 
the chosen time window. Time is measured in reduced units 

m. 
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FIG. 2: Top panel: Density correlators predicted by MCT, 
(l>q'^'^{t), for a one-component SW system for different tem- 
peratures, crossing the glass-glass transition at </) ~ 0.54, for 
qa = 14.5. Such g- value has been chosen to reproduce compa- 
rable differences in the non-ergodicity parameters. The tem- 
peratures have been chosen such that T/Tc {Tc is the MCT 
transition temperature) is close the simulation one for the 
reported T's. Bottom Panel: Same curves from the sim- 
ulation for qaB — 25. Temperatures from top to bottom: 
0.1, 0.2, 0.35, 0.4, 0.5, 0.65, 0.9, 1.5. = 3754. 



+ ty^) = (Pq(*«' + t)pciitw))/{\pciityj)\^) (1) 

where Pq(t) = Si ^^'^ ^ii^) coordi- 

nates of particle i at time t. Fig. ^ shows (f>q{t^ +t,tw) 
for the case T = 0.35, at qasB ^ 25. For > 4000, 
4>qitw+t, tw) becomes independent of for t < 400. We 
confirmed such tu,-independence at all studied T-values. 

Before presenting the numerical data, we discuss the 
ideal-MCT predictions for the one-component system, 
using the Percus-Yevick static structure factor as input, 
extending the calculation reported in Ref. [l7| . Fig. [21 



(top panel) shows 



lMCT 



(t) at different T values. The 



long-time limit fq of the correlation function separates 
two groups of temperatures. At high T values, corre- 
lation functions tend to f^^ « 0.6, while for lower T, 
fq > 0.9, a figure typical of the short-ranged attractive 
glass. The sharp jump in fq is the signature of a discon- 
tinuous glass-glass transition. 

Fig. |21(bottom panel) shows + t^t^) in the i^- 

independent time window, calculated from MD trajecto- 
ries. At the highest T,T = 1.5, (t)q{tw+t, tw) quickly ap- 
proaches a plateau, consistent with the hard-sphere glass 
value, as expected since the kinetic energy is larger than 
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the depth of the square well. On cooling, the approach 
to the plateau takes longer and longer, consistently with 
the MCT predictions. Around T — 0.4, the decay is 
well represented by a logarithmic law, in agreement with 
the location of the MCT higher order singularity — the 
source of the logarithmic dynamics For T < 0.4, 
the correlation functions show a very ill-defined plateau, 
around 0.9, a value consistent with the MCT prediction 
for the short-ranged attractive glass. While in the the- 
oretical predictions the attractive plateau lasts for ever, 
in the simulated systems all correlation functions, inde- 
pendently from T, appear to decay well below 0.9. The 
correlation functions for the two extreme temperatures in 
FigEl (bottom panel) show that a clear difference in the 
dynamics is indeed observed. However, the sharp glass- 
glass transition predicted by ideal MCT is not found. 
These results suggest that the range of stability of the 
attractive glass is considerably shorter than theoretically 
expected. Even lowering T down to 0.1, i.e. when T is 
one tenth of the potential well depth and the system is 
expected to be deep in its attractive glass phase, the at- 
tractive plateau lasts less than two decades in time[27|. 
We also note that the time window in which the correla- 
tion function is close to the attractive plateau value does 
not extend on increasing t^, as clearly shown in FigQ 
The absence of a clear short-ranged attractive plateau 
seems to suggest that the attractive cages are not as sta- 
ble (in time) as the hard-sphere ones and that particles 
do manage to break the bond confinement . 

A possible explanation of the observed instability of 
the attractive glass can be formulated in term of the so- 
called activated (or hopping) processes, decay processes 
which are not included in the ideal MCTj23. While in 
hard-spheres hopping phenomena are not significant due 
to the absence of an energy scale, in molecular systems 
it is well known that the "MCT glass" is destabilized by 
activated processes, which slowly restore ergodicity. In- 
deed, in these systems the calorimetric glass transition 
temperature is located well below the MCT critical tem- 
perature. In short-ranged attractive colloids, activated 
processes can be associated with thermal fluctuations of 
order mq, which are able to break the bonds. These pro- 
cesses generate a finite bond lifetime and destabilize the 
attractive glass. If is sufficiently large, the underly- 
ing hard-core repulsive glass should emerge as limiting 
disordered arrested structure. 

It is tempting to formulate the hypothesis that an ideal 
MCT short-ranged attractive glass would be stable if 
bond-breaking processes were negligible. Such ideal at- 
tractive glass should be characterized by bond-cages with 
infinite lifetime, in analogy with the hard-sphere glass 
which is characterized by neighbor-cages of infinite life- 
time. To substantiate this claim, we artificially build a 
model in which bonds are never broken on the time scale 
of our numerical observation. To do so we add a finite 
barrier of infinitesimal width in the potential, just outside 
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FIG. 3: Comparison of the density correlators (j)q{tw +t,tw) 
for the square well (SW) model and the SW model comple- 
mented by a barrier of infinitesiman width (SWB) and height 
IOOmo, as sketched in the figure, ~ 3754. The dashed lines 
are fits to the SW data, reported as a guide to the eye to help 
comparing the SW and SWB fq values. 



the attractive well (see Fig|2} , with the purpose to make 
the bonds longer lived, and solve the Newton equations 
exactly js^. The nice feature of this modified system is 
that the thermodynamics is identical to the one of the 
original potential, since the width of the repulsive bar- 
rier is infinitesimal. Dynamics are instead different, since 
bond-breaking takes now place on longer time-scales; for 
example, by choosing as height of the repulsive barrier 
IOOuq, we find that no bond is formed or broken in the 
timescale of the simulation. 

Fig. 131 contrasts the decay of ^^(tiu + t,tw) for dif- 
ferent T values in the range of T where attractive glass 
is expected, both for the SW model and for the case 
of the square well complemented by a barrier (SWB) of 
height IOOmq. The same initial configurations (at tw) are 
used for both models, so that at time t = the bonding 
pattern is identical. By construction, only decorrelation 
processes which do not modify the bonding pattern, i.e. 
the rattling motion within the bond-cages, are possible 
in the SWB model. We observe that, in the case of the 
SWB model, the correlation functions decay to a stable 
plateau, whose value, larger than 0.9, is consistent with 
the expected value for the attractive glass. In contrast, 
correlation functions for the unmodified SW model, al- 
ready at short times, decay below the plateau value. Still, 
the extrapolated amplitude of the long-time decay (see 
Fig. |2Jl provides an estimate of fq very similar to the 
one of the SW model plus barrier. These results sup- 
port the view that ideal-MCT, by neglecting the bond- 
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breaking processes, predicts a stability window for the 
short-ranged attractive glass larger than numerically ob- 
served. 

The present study suggests that the ideal attractive 
glass line in short-ranged attractive colloids has to be 
considered, in full analogy with what has been found 
in the study of glass-forming molecular liauids[29j. as a 
cross-over line between a region where ideal-MCT predic- 
tions are extremely good (in agreement with the previous 
calculation in the fluid phase) and an activated-dynamics 
region, where ideal-MCT predictions apply in a limited 
time window. The anomalous dynamics which stem from 
the presence of an higher-order singularity in the MCT- 
eauationsflll still affect the dynamical processes in 
the fluid and in the glass, even if activated processes 
pre-empt the possibility of fully observing the glass-glass 
transition phenomenon, at least in the SW case. Short- 
ranged inter-particle potentials which stabilize bonding 
could produce dynamics which are less affected by hop- 
ping processes, favoring the observation of the glass-glass 
phenomenon^ • 

We acknowledge support from MlUR COFIN 2002 and 
FIRB. We thank S. Buldyrev for the MD code, W. Gotze 
for a critical reading of the manuscript and M. Fuchs and 
A.M. Puertas for discussions. We acknowledge support 
from INFM Iniziativa Calcolo Parallelo. 



[1] P. N. Pusey, A. D. Pirie and W. C. K. Poon, Physica 
A 201, 322 (1993); S. M. Ilett et al, Phys. Rev. E 51, 
1344 (1995); P. N. Segre et al, Phys. Rev. Lett. 86, 6042 

(2001) ; V. Trappe et al, Nature 411, 772 (2001). 

[2] H. N. W. Lekkerkerker et al, Europhys. Lett. 20, 559 
(1992); P. Bolhuis, M. Hagen abd D. Frenkel, Phys. Rev. 
E 50, 4880 (1994); C. F. Tejero et al, Phys. Rev. E 51, 
558 (1995); V. J. Anderson and H. N. W. Lekkerkerker, 
Nature 416, 811 (2002); K. A. Dawson, Curr. Opin. Col- 
loid Interf. 7, 218 (2002). 

[3] F. Sciortino, Nature Materials 1, 145 (2002). 

[4] In experiments, depletion interactions control the in- 
terparticle potential. The role of T is provided by the 
amount of non-adsorbing polymer, while the range of in- 
teraction is controlled by the length of the polymer. 

[5] F. Mallamace et al, Phys. Rev. Lett. 84, 5431 (2000). 

[6] K. N. Pham et al, Science 296, 104 (2002). 

[7] T. Eckert and E. Bartsch, Phys. Rev. Lett. 89, 125701 

(2002) . 

[8] W. R. Chen, S. H. Chen and F. Mallamace, Phys. Rev. 

E 66, 021403 (2002). S. H. Chen, W. R. Chen and F. 

Mallamace, Science, 300, 619 (2003). 
[9] G. Foffi et al, Phys. Rev. E 65, 050802(R) (2002). 
[10] E. ZaccareUi et al, Phys. Rev. E 66, 041402 (2002). 
[11] F. Sciortino, P. Tartaglia, E. ZaccareUi, 

cond- mat/0304192 (2003). 
[12] A. M. Puertas, M. Fuchs and M. E. Cates, Phys. Rev. 

Lett. 88, 098301 (2002); Phys. Rev. E 67, 031406 (2003). 
[13] W. Gotze in Liquids, Freezing and Glass Transition 



edited by J. P. Hansen , D. Levesque D, and J. Zinn- 

Justin (Amsterdan: North-Holland) p 287, 1991. 
[14] W. van Megen and S. M. Underwood, Phys. Rev. Lett. 

70, 2766 (1993); Phys. Rev. E 49, 4206 (1994). 
[15] L. Fabbian et al, Phys. Rev. E R1347 (1999), and Phys. 

Rev. E 60, 2430 (1999). 
[16] J. Bergenhohz and M. Fuchs, Phys. Rev. E 59, 5706 

(1999). 

[17] K. A. Dawson et al, Phys. Rev. E 63, 011401 (2001). 
[18] E. ZaccareUi et al, Phys. Rev. E 63, 031501 (2001). 
[19] W. Gotze and M. Sperl, Phys. Rev. E 66, 011405 (2002). 
[20] W. Gotze and M. Sperl, J. Phys. Condens. Matt. 15, 
S869 (2003). 

[21] Lengths are measured in units of as, time in units of 
(JB{Tn/u(,Y^^ and fcs, the Boltzmann constant, is fixed 
to 1. Energy and T are thus measured in units of uo. The 
packing fraction (j) is defined as 7r/6(A^AO"^ + Nbo-%)/V , 
where Na = Nb = 350 and V is the volume. 

[22] M. Sperl, preprint (2003). 

[23] It is known that MCT overestimates the tendency to form 
a glass. A mapping is required to compare theoretical and 
numerical results. In the present bilinear mapping 

i4> 1.897 4> - 0.3922, T 0.5882 T - 0.225) success- 
fully superimposes the MCT predictions for the two ideal 
glass lines onto the corresponding numerical estimates, as 
discussed at length in R,ef.|lll|. 

[24] To quickly compress, we run a simulation of a system of 
particles interacting with a hard-core followed by a re- 
pulsive shoulder. The configuration to be compressed is 
used as initial conditions in the Newton Equations. The 
heigth of the shoulder is set to a value much larger than 
T, forcing closeby particles to quickly separate. As soon 
as the energy of the system goes back to zero, the sim- 
ulation is interrupted and the configuration is saved. A 
proper rescaling of the coordinates produce a configura- 
tion consistent with the original square- well potential but 
with a larger packing fraction (fixed by the width of the 
repulsive shoulder) . In this way, we were able to generate 
configurations up to about (j> — 0.65. 

[25] It is worth observing that after a waiting time of quench 
of about 250 MD units, no significant release/absorption 
of heat is observed in the studied time window. This is 
due to the high packing fraction and the associated slow- 
ness of the particle dynamics which allows very minute 
configurational changes. 

[26] W. Kob and J.-L. Barrat, Phys. Rev. Lett. 78, 4581 
(1997); L. F. CugUandolo, condmat/0210312. 

[27] The presence of a short-lifetime of the attractive glass 
cage was also encountered in p^.ll^. 

[28] Despite the fact that the lifetime of the attractive glass 
is not infinite, the presence of a short-ranged bonding 
significantly affects the time evolution of (pq and, corre- 
spondingly, the frequency dependence of the elastic mod- 
ulus on varying T. 

[29] S. P. Das and G. F. Mazenko, Phys. Rev. A 34, 2265 

(1986) ; W. Gotze and L. Sjogren, Z. Phys. B 65, 415 

(1987) . 

[30] We have implemented an event-driven algorithm, where 
in addition to hard-core collisions, the events include also 
the barrier crossing as explained in D. C. Rapaport, The 
Art of Molecular Dynamic Simulation, Cambridge Uni- 
versity Press, 1995. 



